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ABSTRACT 

Wc present a H I 21 cm absorption survey with the Green Bank Telescope (GBT) of galaxy-quasar 
pairs selected by combining galaxy data from the Sloan Digital Sky Survey (SDSS) and radio sources 
from the Faint Images of the Radio Sky at Twenty- Centimeters (FIRST) survey. Our sample consists 
of 23 sightlincs through 15 low-redshift foreground galaxy - background quasar pairs with impact 
parameters ranging from 1.7 kpc up to 86.7 kpc. We detected one absorber in the GBT survey 
from the foreground dwarf galaxy, GQ1042+0747, at an impact parameter of 1.7 kpc and another 
possible absorber in our follow-up Very Large Array (VLA) imaging of the nearby foreground galaxy 
UGC 7408. Both of the absorbers are narrow (FWHM of 3.6 and 4.8 km s _1 ), have sub Damped 
Lyman a column densities, and most likely originate in the disk gas of the foreground galaxies. We 
also detected H I emission from three foreground galaxies including UGC 7408. Although our sample 
contains both blue and red galaxies, the two H I absorbers as well as the H I emissions are associated 
with blue galaxies. We discuss the physical conditions in the 21 cm absorbers and some drawbacks of 
the large GBT beam for this type of survey. 
Subject headings: galaxies: abundances — galaxies: ISM — quasars: absorption lines 



1. INTRODUCTION 

Inflow and outflow of gas between the intergalactic 
medium (IGM) and the disks of galaxies is a crucial as- 
pect of galaxy evolution. In the last decade our theo- 
retical understanding on how galaxies acquire gas has 
advanced considerably. Traditionally, gal axies were be- 
lieved to accrete via the "hot" mode (jWhite fe Reesl 
Il978t IWhite fc FrenldfT991| ). where accreting gas shock- 
heats to the virial temperature and then all halo gas 
within a "cooling radius" cools to form stars. Recent 
hydrodynamic cosmological simulations have discovered 
that gas accretion also occurs in the so called "cold" 
mode where the gas is able to cool as it descends into 
galactic potentials and thus never approaches the virial 
temperature (see iBirnboim fc Dekel 120031 : iKeres et al.l 
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Foundation operated under cooperative agreement by Associated 
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scope, which is owned and operated by the Astrophysical Research 
Consortium. 



l2005f ). Unlike the hot-mode, the cold- mode accretion 
occurs through long filaments connecting the extended 
disks and halos of galaxies to their surrounding in- 
tergalactic medium (IGM). Although initial results by 
iKeres et al.l (|2005() snowed that the cold mode accretion 
is dominant in galaxies with lower masses ( < 10 103 M^) 
at red shift z = 0, recent simulations by iKeres et al.l 
( 2009) have shown that cold-mode accretion also persists 
even in massive galaxies with hot envelopes. In addition 
to filamentary accretion of cool gas from the intergalac- 
tic medium, cool gas is al so expected to co n dense out 
of the hot galactic halos. iMaller fc Bullockl (|2004ft ar- 
gued that small clouds of cool gas condense out of hot 
gaseous halos through a multiphase process due to ther- 
mal instabilities. These halo clouds are predicted to have 
a characteristic size of « 1 kpc and are expected to be 
another continuous source of cool gas for galaxies. 

On the other hand, a substantial amount of material 
is ejected from galactic d isks into their halos via galactic 
winds (see the review bv lVeilleux et al.ll2005h . Evidence 
of outflows has been detecte d in low-ionization metal 
absorption lines of Na I (e.g. ISchwartz fc Martini l2004t 
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Rupke et al I l2005allH) and Mg II jlTremonti et all 120071; 






Rubin et al.ll2010t ). iRubin et all fl2010|) found that the 
strength of the outflow correlates with the absolute star- 
formation rate (SFR) and claimed that the mass outflow 
rates were of the order of the SFRs for the i r samp le. As- 
suming a negligible halo drag. lRupke et ail (|2005b| ) found 
that only 5-10% of the neutral material in starburst- 
driven winds escape into the IGM. In other words, the 
wind speeds are not larger than the galactic escape veloc- 
ity for most galaxies and thus the outflowing material is 
primarily retained in the hal os of galaxies. Accordi ng to 
the Galactic Fountain model (jShapiro fc Fieldll"l976l ). the 
baryonic material deposited by the winds can eventually 
cool down and rain bac k into the galax ies. This model 
and subsequent work bv lBregmanl (|1980h predicts the ex- 
istence of high-velocity clouds (HVCs) as a by-product 
of the cooling of the hot gas. 

Besides inflow and outflow signatures, there is growing 
evidence of the extended nature of the galactic disks. 
Mg II absorbers, which trace cool low-ionization gas, 
have been believed to be associated with extended disks 
(lWagonedll967l:lBowen et~all[l995l: ICharlton fc ChuTchil 
Il996| ). Observations of sightlines probing foreground 
galaxies at impact parameters up to 75 h~ l kpc with in- 
clination angles between 40° to 75° have co nfirmed Mg II 
absorbers exhibiting disk- li ke rotation (jSteidel et al.l 
[2001 IKacprzaket"all[20T0l) . ISteidel et al.l ((2002T) found 
that four of five background QSOs that lie to one side 
of a foreground galaxy exhibit Mg II absorption from 
the latter that are shifted in the same sense from the 
systemic redshift of the foreground galaxy. They con- 
cluded that the kinematics of their sample could be un- 
derstood in terms of a disk-like rotation combined with 
a lagging halo. This w as confirmed in a recent study by 
iKacprzak et al.l (|2010D . where most of the Mg II absorp- 
tion systems in their sample fully reside to one side of the 
galaxy systemic velocity. They also found a close cou- 
pling between galaxy inclination and Mg II absorption 
velocity spread. While it is clear that there is cool gas at 
these distances ( > 100 kpc for h = 0.7 km s _1 Mpc -1 ), 
it is unclear how much baryonic mass exists in these sys- 
tems. Part of the problem is that these observations do 
not probe the hydrogen directly, and the coupling of neu- 
tral gas and metals is not well understood. 

To properly characterize the total baryonic content in 
the extended disks as well as gaseous halos, H I ob- 
servations must be made at a wide range of distances 
from the galaxies. Absorption spectroscopy is a pow- 
erful tool for probing low-column density systems with 
high sensitivity independent of the redshift of the ab- 
sorbing material. Unfortunately, the Lyman a tran- 
sition for z < 1.7 galaxies occurs in the ultraviolet 
(UV) band and thus ground-based observations of these 
systems are not possible. However, the H I 21 cm 
hyperfine transition for low— z systems can be ob- 
serve d with existing radio facilit i es. Various authors 
(e.g.. iKanekar fc Chengaiurlll997t lLane fc Briggl 120011: 



Vermeulen et al.ll2003t iDarling et al.l l2004t iKeenev et all 



20051: iGupta et al.H2007L 12003 . l2010h have probed neu- 
tral gas in extragalactic sources at various redshifts us- 
ing 21 cm absorption. I n their recent study of the 
radio galaxy Centaurus A, IStruve et al.l (|2010f ) discov- 
ered an H I 21 cm absorbing cloud of column density 



~ 10 20 cm -2 (assuming T s = 100 K and filling factor 
1) in a region beyond the optical extent of the galaxy, 
at a distance of 5.2 kpc from the galactic nucleus us- 
ing the northeastern radio lobe of the galaxy itself as a 
background source. While there is a substantial amount 
of literature on the using H I absorption toward ra- 
dio lobes to probe the immediate vicinity of AGN (e.g. 
ISaikia fc Jamro zv 2009]), this example shows how radio 
galaxies can also be used to probe the ISM of galaxies 
far from the nuclear region." 

Another advantage of radio absorption studies is that 
radio sources are often extended and thus unlike UV or 
optical absorption spectroscopy, which is a single pen- 
cil beam study, absorption in the radio band provides a 
unique opportunity to probe multiple sight lines through 
the ISM of foreground galaxies simultaneously, and thus 
is an excellent way to map the small-scale structures in 
extragalactic ISMs. Such small-scale structures down to 
AU-scales were detected in th e Milky Way as early as 
in 1976 bv iDieter et al.l (J1976I ) and later confirmed by 
iDiamond et al.l ( 19891) . Over the last few years, multi- 
ep och mapping of the e xten ded radio emission of quasars 
bv lBrogan et al.1 (|2005[ ) and lLazio et~aT1 (|2009l ) have con- 
firmed the small-scale structures in the ISM of the Milky 
Way detected via 21 cm absorption. These authors re- 
port that the covering fraction of these structures are 
roughly 10% which implies that their volume filling frac- 
tion is only about 1%. Similar spatial variations in the 
ISM of galaxies outside the Milky Way have not yet been 
detected. One way to constrain the sizes and physi- 
cal nature of the cold clouds, including their spin tem- 
perature and covering fraction, is to conduct a survey 
of H I absorption against multiple bright background 
sources that probe a single galaxy. Such a survey would 
provide unique information on the (1) structure of the 
cold gas clouds in a galaxy, (2) variation of cloud prop- 
erties as a function of distance from the galaxy/galactic 
disk; and (3) the cold gas distribution as a function of 
global properties (such as color, stellar mass) of galaxies. 

In this paper, we investigate the existence and prop- 
erties of cold gas in the disks and halos of low-redshift 
galaxies (z < 0.3). We seek to constrain the frequency 
and strength of the H I distribution as function of the 
distance from the center of a galax y. A similar study was 
carried out by IGupta et al.l (|2010t ) with 5 quasar-galaxy 
pairs using the Giant Metre Radio Telescope (GMRT). 
The most important difference between their observa- 
tional set up and ours is in the spectral and spatial res- 
olution. The two studies are highly complementar y. On 
one hand, the study conducted by IGupta et al.l has a 
high spatial resolution which enables them to detect ab- 
sorption more efficiently. On the other hand, our study 
achieved a high spectral resolution, which is crucial in 
identifying narrow H I features. The 21 cm absorb- 
ing gas is expected to be cold and thus the line widths 
can be extremely narrow. Such fe atures were observed 
in th e Milky Way (~1 km s _1 by iKnapp fc Verschuurl 
I972T) . the Large M agellanic Cloud (LMC) (1.4 km s" 1 by 



Dickev et al.lll99l. and Small Ma gellanic Cloud (SMC) 
(1.2 km s _1 bv lDickev et al.|l2000D Unlike the extended 



sample used bv IGupta et al.l (|201C0 in their analysis, we 
use an unbiased sample, which was selected solely on the 
basis of the impact parameter without any consideration 
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to absorption detection through other line transitions. 

We also wish to investigate the H I distribution in 
terms of galaxy luminosity, mass, and color. This goal 
is motivated by the fact that the intensity of accretion 
and outflow processes are expected to vary from galaxy 
to galaxy. For instance, blue star-forming galaxies are 
expected to have a significant reserve of gas compared 
to red galaxies. Similarly, cold-mode accretion is pre- 
dicted to be predominant in galaxies of baryonic mass 
< 10 103 Mq. Secondly, we intend to study the proper- 
ties of H I clouds, such as their size, filling factor, and 
temperature, in various regions of an individual galaxy. 
We begin by presenting our sample and elaborating on 
the selection procedure, followed by our HI 21 cm ob- 
servations in section 2. In section 3 we present our re- 
sults and discuss their implications in section 4. Finally, 
we summarize our findings and conclusions in section 5. 
Throughout the paper we use H = 70 km s _1 Mpc -1 , 
Q m = 0.3, and ft A = 0.7. 

2. OBSERVATIONS 

2.1. Sample Selection 

Our sample was selected in a two-stage procedure. 
First, a set of background QSO sightlines close to 
foreground galaxies were identified by cross-correlating 
galaxy positions with quasar positions from objects ob- 
served by the Sloan Digital Sky Survey (SDSS). The 
SDSS has redshift information for some of the galaxies, 
while others have photometric redshift determinations. 
Hence to identify quasar-galaxy pairs at various impact 
parameters from the entire SDSS catalog, we used one of 
the four criteria listed below. 

1. For SDSS galaxies with spectroscopic redshifts be- 
tween 0.005 < Zg a i < 0.2, we searched for SDSS 
QSOs that had sightlines passing within 20" of 
a galaxy, or, within 1.2 x rp(r), where rp(r) is 
the r— band Pctrosian radius of a galaxy in arcsec. 
This latter selection was useful for finding low red- 
shift galaxies where the QSO-galaxy impact param- 
eter was > 20" but where the sightline still passed 
close to the optical extent of a galaxy. Similarly, 
we also searched for QSO sightlines from galax- 
ies separated by < 30" but where the galaxy had 
rp(r) > 8.0, in order to find very low redshift galax- 
ies (with, therefore, large rp(r) on the sky) and 
where larger angular separations still translated to 
relatively small physical impact parameters. In or- 
der to remove galaxies whose absorption features 
might be confused with the AGN activity of the 
QSO, we excluded galaxies whose redshifts were 
within Az = 0.1 of the QSO's redshift. 

2. We searched for SDSS QSO sightlines that passed 
within 30" of SDSS galaxies that had no measured 
redshifts, with the intention of obtaining redshift 
information using available telescopes. We discuss 
these data in Section 12.21 



4. For probing the smallest impact parameters, we 
chose foreground galaxies that were identified by 
detecting emission lines su perimposed on the back- 



ground QSO spectru m (|Quashnock et al.l 120081 : 
iBorthakur et al 
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3. We also searched for SDSS QSOs within 1.5 D 25 
of all galaxies in the Third Reference Catalogue 
of B right Galaxies (RC3) (|de Vaucouleurs et al.1 
119911 ). where D25 is the major isophotal diameter 
of a galaxy measured to a surface brightness level 
of [iB = 25.0 square mags. 



Based on the above criteria, we selected about a thou- 
sand QSO-galaxy pairs. In the second stage of selecting 
suitable QSO-galaxy pairs, we extracted the radio flux 
densities associated with the optical QSOs from the Faint 
Images of the Radio Sky at Twenty- Centimeters (FIRST) 
Survey ([Becker et al.lll995T ). We chose a limiting 1.4 GHz 
flux density value of 25 mJy for the background quasar, 
corresponding to 3 a optical depth of 0.35 for 2 hours of 
on-sourcc integration for our observational settings de- 
scribed in the next section. Based on the observable red- 
shif t range and exclud ing pairs that were being observed 
by ([Gupta et al.ll2010D with the GMRT, this flux density 
cut resulted in a sample of 18 galaxy-quasar pairs. 

We carried out observations of the 18 galaxy-pairs 
with the Green Bank Telescope (GBT). We encountered 
severe radio frequency interference (RFI) for three of 
the pairs, thus reducing our sample to 15 galaxy-quasar 
pairs. Some of the background quasars are extended ra- 
dio sources with multiple peaks at the resolution of the 
FIRST Survey, and therefore allow us to probe different 
regions of the foreground galaxies. The GBT spectrum 
for such a pair would have contributions from each of the 
peaks seen in FIRST image as the GBT beam encom- 
passes all the peaks of the extended quasars. Although, 
this would cause confusion in the case of a detection, 
however, this is particularly advantageous for nondetec- 
tions. More stringent limits can be estimated for each 
of the background continuum peaks for nondetections. 6 
of our 15 galaxy-quasar pairs have extended quasars and 
these GBT nondetections essentially probe 14 slighlincs. 
Hence, we present our data for the 23 sightlines through 
15 galaxy-quasar pairs. 

Details on the sightlines along with their 1.4 GHz ra- 
dio flux densities are presented in Table [TJ Radio con- 
tours from the FIRST survey overlaid on the composite 
SDSS r, g, and i-band images for our sample are pre- 
sented in Figure [TJ The QSO is centered in each im- 
age and the foreground galaxy is marked by a white ar- 
row. The QSO name is printed at the top left-hand side 
of the image whereas the sightline numbers used in Ta- 
blcUJare reproduced near the galaxy or background radio 
source. Including multiple radio peaks, our sample covers 
a wide range of impact parameters from 1.7 to 86.7 kpc as 
well as a wide variety of galaxy properties (see next sec- 
tion). The sightline toward the background QSO SDSS 
J104257. 58+074850. 5 allows us to probe the inner stellar 
disk of the foreground dwarf spiral gala xy GQ1042+0747 
at an impact parameter of 1.7 kpc (see IBorthakur et al.l 
(|2010t ) for details). Besides a wide continuous range of 
impact parameters, we also cover a wide range of galaxy 
properties including luminosity, mass, and color. 

2.2. Optical Properties and Galaxy redshifts 

For foreground galaxies without redshift information 
from SDSS, we conducted spectroscopic observations 
from the Apache Point Observatory (APO) to obtain 
their redshifts. The APO redshifts were measured ei- 
ther by fitting Gaussian profiles to emission lines or from 



TABLE 1 
A Detailed Description of Our Sample with Information on the Background Optical QSOs, Their Radio Counterparts, and the Foreground Galaxies. 
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1 


J010644.15-103410.5 


010644.124-103410.55 


265.90 


J010643. 94-103419. 3 


0.1460 


APO 


18.68 
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10.0 


Blue 


23.4 


2 
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J074841. 77+173456.6 
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102258.415+123426.26 


93.98 
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APO 
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SDSS 
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125247.588+474042.81 
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1.861 


11.1 
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45.06 
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APO 
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11.1 
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0.0478 


SDSS 
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SDSS 
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20 
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70.21 
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SDSS 
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21 
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J115839.90+625428.0 


115842.780+625442.26 


69.64 


J115840.67+625423.7 


0.2596 


SDSS 


18.49 


2.887 


11.1 


Blue 


94.5 


a Q 


upta et al. 2010 in press 
cntificd through emission 


lines on the QSO spectri 




















h Id 


nil (seelQuashnock et al.H2008l; IBorthakur et al.H2010h 





Also known as UGC 7408. 



-t 



GBT Absorption Survey for H I in Halos of Galaxies 



5 o 




5 



15 

RA Offset " 




15 -15 -30 

RA Offset " 






30 


J0910+4636 




Mk 




15 




1 


* 


3 







1 


1 




-15 
-30 




7 


•m 



15 -15 -30 

RA Offset " 



15 

RA Offset ' 




15 -15 -30 

RA Offset " 





15 -15 -30 

RA Offset " 



15 

RA Offset " 



90 60 30 -30 -60 -90 
RA Offset " 




15 

RA Offset " 




J1328+6221 






19 


I 



15 

RA Offset " 



15 

RA Offset ' 



Fig. 1. — SDSS composite image of i-, r-, and g-band images for our sample of 15 galaxy-quasar pairs. The quasar is at the center of the 
image and its name is printed at the top left side of each plot. The yellow contours shows the 1.4 GHz VLA B-connguration continuum 
emission from the FIRST survey at levels 1, 5, 25, and 125 mjy/beam. The foreground galaxies arc marked with a white arrow and each of 
the sightlines is numbered. Due to the large extent of the background quasars in 1.4 GHz some of our systems covered multiple sightlines. 
Details of our sample are presented in table [T] 
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Fig. 1. — Continued. 



TABLE 2 

apo redshifts for foreground galaxies using emission 

and Absorption lines. 
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z a 

abs — line 


5z a , , . 

a bs — in i a 


J010643. 94-103419. 3 


0.1460 


- 


- 


J084913. 80+275735.9 


0.2356 


0.2354 


0.0004 


J084912. 42+275740.4 


0.1941 


0.1948 


0.0008 


J091010.55+463633.3 


0.0998 


0.0998 


0.0004 


J102257.92+123439.1 


0.1254 


0.1253 


0.0004 


J110736.88+090113.6 


- 


0.1050 


0.0003 


J132839. 89+622136.0 


0.0420 


0.0423 


0.0005 


J141629. 25+372120.4 


0.0341 


- 


- 



a Estimated using IRAF fxcor routine. 

the stellar absorption line measurements made using the 
Image Reduction and Analysis Facility (IRAF) "fxcor" 
routine, and cross-correlating the blue spectra with that 
of a radial velocity standard. The redshifts are presented 
in the sixth column of Table Q] with more details in Ta- 
ble El With the exception of two galaxies, the remain- 
ing APO redshifts arc measured using absorption-lines. 
The r-band magnitude and luminosity for the foreground 
galaxies are presented in the eighth and ninth columns 
of Table [TJ The luminosities, in units o f L* , were es- 
timat ed using M; = -21.2 (Table 2 of iBlanton et al.l 
120031 ) and were k-corrected for r-band absolute magni- 
tude at z = 0.0 . The k-cor rections were derived us - 
ing the k-correct algorithm by IBlanton fc Roweisl (|2008f ) 
that uses SDSS photometric u, g, r, i, and z-band data. 
Due to the unavailability of SDSS photometric data for 
GQ1042+0747 we have used g, r, and i— band photomet- 
ric values from lBorthakur et al.l (|2010ft to estimate the k- 
correction for this galaxy. We also present stellar masses 
for the foreground galaxies, which were e stimated using 
Equation (1) from [Mcinto sh et al.l (l2008f ) based on the 
k-corrected magnitudes. 

From the galaxy color-magnitude diagram that re- 
lates galaxy luminosity and color one can classify galax- 
ies as part of the red sequence or the blue cloud. 
For our sample, we follow the empirical results from 
Ivan den Bosch et al.l (J2008L equation A6), which uses the 
parameter °-°A defined as 

°-°A = ao (g -r)- 0.10 (log[M„] - 10.0) (1) 

where 00 (g — r) is based on g and rband magnitudes that 
have been k-corrected to z = 0.0 and M* is the stellar 
mass of the galaxy. The two populations are found to 
have an overlap in terms of the 00 A value intersecting 



at 00 A = 0.675 (see Figure Al of Ivan den Bosch et al.l 
2008). Hence we use this value to distinguish the two 
populations and we classify galaxies with 0,0 A < 0.675 
as belonging to the blue cloud and ones with 00 A > 0.675 
as red sequence. In total, our sample consists of nine blue 
cloud and six red sequence galaxies. 

2.3. Radio Observations 

In order to search for cold gas in the outer disks and 
halos of galaxies, we carried out GBT observations of our 
sample followed by higher spatial resolution Very Large 
Array (VLA) and Very Large Baseline Array (VLBA) ob- 
servations. Details of the VLA and VLBA observations 
of the pair SDSS J10 4257.58+074850.5 - GQ10 42+0747 
have been presented in iBorthakur et al.l (|2010f ) and are 
not presented in this section. 

2.3.1. 21 cm Spectroscopy with the GBT 

We conducted GBT observations of the HI 21 cm 
and OH 18 cm transitions for 18 galaxy-quasar pairs 
(Tabled]) under project IDs GBT06B-052, GBT08A-082, 
and GBT09C-039 in August 2006, March through July 
2008, and October and November 2009 respectively. We 
used the dual polarization L-band system with 9-level 
sampling. The intermediate frequency (IF) was set to 
yield a channel width of 1.56 kHz (0.33 km s -1 ) using 
8192 channels over a total bandwidth of 12.5 MHz. Ob- 
servations were made in the standard position-switching 
scheme by cycling through the ON-OFF sequence, then 
dwelling for 300 s at each position. Each 300 s scan con- 
sisted of thirty 10 s integrations to minimize data corrup- 
tion by sporadic RFI. The OFF position was chosen +20' 
offset in Right Ascension from each of the sightlines. This 
was done to track the possible presence of a source in the 
OFF position. During each of the observing sessions, we 
used one of the three standard flux calibrators - 3C 48 
(16.5 Jy), 3C 147 (22.5 Jy), and 3C 286 (15.0 Jy) - for 
pointing and estimating antenna gain. Local pointing 
corrections (LPC) were performed using the observing 
procedure AutoPcak and the corrections were then au- 
tomatically applied to the data. This should result in 
a pointing accuracy of 3". Similarly, flux calibrations 
were performed by applying antenna gain to data for 
each session separately (values are presented in Table [3]) . 
The corrections also include attenuation due to opacity 
of the air (around 1%). 
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The data were analyzed using the NRAO software 
GBTIDL. Most of the sightlines were observed for al- 
most 2 hours on-source. However, due to the corruption 
of the data by RFI the noise levels achieved for each of 
the sightlines span a wide range. Severe RFI corrupted 
entire datasets for sightlines toward the background radio 
QSOs 082153+503120 at z = 0.1640, 084914+275729 at 
z=0.2357, and 115836+625415 at z=0.2596 rendering the 
data unusable. Besides RFI, the GBT data also suffered 
from standing waves for some of the brighter sources, 
which limited our H I emission detection capabilities, 
especially for H I features with A V > 300 km s _1 . 
The details of the noise properties and GBT measure- 
ments are presented in Table [3] No OH emission was 
detected in any of our sources and rms noise levels are 
presented in the above mentioned table. 

2.3.2. Very Large Array Observations of J 122115+454838 

In order to understand the effect of the large GBT 
beam size on the survey, we imaged the H I distribution 
in the foreground galaxy J122115+454838 (UGC 7408) 
with the VLA under program AY190 on 28 June and 
16 July 2008. The observations were carried out in the 
VLA D-configuration for 3.5 hours of on-source integra- 
tion. We used dual polarization using 2 intermediate 
frequency (IF) correlator modes to achieve a total band- 
width of 1.5625 MHz (330 km s" 1 ) with a 12.2 kHz res- 
olution after Hanning smoothing. This correlator setup 
allowed us to cover the entire H I line-width with a ve- 
locity resolution of 2.6 km s _1 . The data were cali- 
brated following the standard VLA calibration and imag- 
ing procedures using NRAO's Astronomical Image Pro- 
cessing Software (AIPS). The data were imaged using 
both natural and uniform weighting schemes. The nat- 
ural weighting was applied to enhance the fainter emis- 
sion features whereas uniform weighting was applied to 
achieve a higher spatial resolution for the purpose of ex- 
tracting a spectrum with minimum contamination from 
the H I emitting region. The synthesized beam pro- 
duced using natural weighting was 67.51" x 45.31" (Po- 
sition angle= —86.75°) and uniform weighting was 52.17" 
x 38.11" (Position anglc= -80.56°). The absolute un- 
certainty in the resulting flux density scaling is about 
10%, and this is the formal uncertainty we quote for all 
physical parameters derived from the flux density. The 
average rms noise achieved, after combining the two po- 
larizations, is 1.07 mjy /beam/channel. 

3. RESULTS 

We detected H I emission from three fore- 
ground galaxies: J122115.22+454843.2 (UGC 7408), 
J141629.25+372120.4, and J160659. 13+271642.6, thus 
suggesting the presence of gas-rich galaxies in our sample. 
Figure [2] shows the smoothed GBT H I spectra for these 
galaxies. The spectra are presented in terms of the op- 
tical velocity of the galaxies, which show very good cor- 
respondence with the H I distribution. The H I masses 
presented in Table [3] were estimated using the relation- 
ship M(ffl) = 2.36xl0 5 J D 2 (S'AU) M Q , where D is the 
luminosity distance to the galaxy in megaparsec (Mpc) 
and the SAV is the velocity integrated H I flux density 
in Jy km s _1 . All the three galaxies are "blue" (as de- 
fined in Section r2.2[) although they vary in luminosity and 
stellar mass. While J122115.22+454843.2 (UGC7408) 
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Fig. 2. — HI spectra for foreground galaxies UGC 7408 (top), 
J141629+372120 (center), and J160659+271642 (bottom). The 
spectra are presented at the rest frame of the optical galaxies to 
show the match between the optical and H I velocity distribution. 
All three spectra have been smoothed to 3 km s — 1 (i.e. 10 pixels). 
The H I mass associated with each of these sources are presented 
in Table [3] 
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TABLE 3 
GBT 21 cm Observations and Measurements of H I Optical Depth and Mass Associated With the Foreground Galaxies. 



SL 



^DvJ ra( ^i 



Si. 4 



Gain 



&HI 



u OH Foreground Galaxy 



Distance 



HI Mass 



1 


010644.124-103410.55 


265.90 


1.68 


3.65 


< 0.02 


2.02 


J010643.94-103419.3 


692.4 


- 


2 


074841.773+173456.82 


41.97 


1.72 


2.45 


< 0.10 


1.63 


J074842. 58+173450.6 


235.2 


- 


3 


074841.786+173512.20 


21.94 


1.72 


2.45 


< 0.19 


1.63 


J074842. 58+173450.6 


235.2 


- 


1 


074842 .084+1 73443. 37 


16.08 


1.72 


2.45 


< 0.26 


1.63 


J074842. 58+173450.6 


235.2 


- 


5 


082153.833+503120.57 


53.94 


1.70 


4.96 


< 0.16 


2.41 


J082153.76+503125.7 


890.5 


- 


6 


084914.282+275729.90 


53.94 


1.72 


6.36 


< 0.20 


1.98 


J084912.42+275740.4 


951.7 


- 


7 


091011.016+463617.87 


163.10 


1.71 


1.31 


< 0.01 


RFI 


J091010.55+463633.3 


459.3 


- 


8 


102258.415+123426.26 


93.98 


1.70 


1.86 


< 0.03 


1.88 


J102257.92+123439.1 


586.4 


- 


9 


102258.552+123439.94 


24.67 


1.70 


1.86 


< 0.13 


1.88 


J102257.92+123439.1 


586.4 


- 


10 


104257.598+074850.60 


381.57 


1.65 


1.2 


0.04962 


RFI 


GQ1042+0747 b 


145.9 


- 


11 


110736.607+090114.72 


36.97 


1.82 


3.81 


< 0.18 


RFI 


J110736.88+090113.6 


484.9 


- 


12 


122105.480+454838.80 


53.69 


1.74 


3.49 


< 0.11 


1.90 


J122115.22+454843.2 C 


8.4 


1.4 


13 


122106.854+454852.16 


21.46 


1.74 


3.49 


< 0.28 


1.90 


J122115.22+454843.2 C 


8.4 


1.4 


11 


122107.811+454908.02 


12.82 


1.74 


3.49 


< 0.47 


1.90 


J122115.22+454843.2 C 


8.4 


1.4 


15 


125247.588+474042.81 


78.72 


1.86 


2.64 


< 0.06 


RFI 


J125249.20+474105.9 


140.9 


- 


16 


125249.326+474042.19 


45.06 


1.86 


2.64 


< 0.10 


RFI 


J125249.20+474105.9 


140.9 


- 


17 


132534.240+495348.47 


26.63 


1.64 


5.90 


< 0.38 


RFI 


J132535. 13+495324.4 


212.2 


- 


18 


132534.565+495342.26 


15.87 


1.64 


5.90 


< 0.64 


RFI 


J132535. 13+495324.4 


212.2 


- 


19 


132840.599+622136.65 


144.10 


1.78 


3.02 


< 0.04 


2.13 


J132839.89+622136.0 


187.0 


- 


20 


134528.765+034720.09 


70.21 


1.91 


3.12 


< 0.08 


RFI 


J134525.30+034823.8 


142.7 


- 


21 


141631.039+372203.01 


30.71 


1.87 


2.75 


< 0.16 


RFI 


J141629.25+372120.4 


149.9 


23.2 


22 


141630.672+372137.09 


30.53 


1.87 


2.75 


< 0.16 


RFI 


J141629.25+372120.4 


149.9 


23.2 


23 


160658.315+271705.86 


141.01 


1.76 


2.34 


< 0.03 


RFI 


J160659. 13+271642.6 


204.8 


81.3 


- 


082153.833+503120.57 


53.94 


RFI 


RFI 


RFI 


1.88 


J082152. 85+503128.4 


815.1 


- 


- 


084914.282+275729.90 


53.95 


RFI 


RFI 


RFI 


2.40 


J084913. 80+275735.9 


1177.1 


- 


- 


115836.882+625415.06 


401.20 


RFI 


RFI 


RFI 


2.20 


J115840.67+625423.7 


1315.4 


- 


- 


115837.852+625420.28 


180.07 


RFI 


RFI 


RFI 


2.20 


J115840.67+625423.7 


1315.4 


- 


- 


115842.780+625442.26 


69.64 


RFI 


RFI 


RFI 


2.20 


J115840.67+625423.7 


1315.4 


- 


a In 


tegrated optical depth binned over 1 km s — 1 bins i.e 3 pixels. The 1 
entified through emission lines on the QSO spectrum (sec Quaslmoc 


imiting values show 3 <7 1 


optical depth. 




b Id 


k et al.ll200St IBorthakur et al.H2010h 





1.0, 



c Also known as UGC 7408 

is a gas rich dwarf galaxy of luminosity 0.003 L* and 
H I mass of 1.4 x 10 8 M , J141629.25+372120.4 
and J160659. 13+271642.6 are ~ 0.1L t galaxies and have 
H I masses of 2.3 x 10 9 M and 8.1 x 10 9 M^ 
respectively. 

The VLA D-configuration image of H I emission asso- 
ciated with UGC 7408 (green contours), overplotted on 
an SDSS r-band image, is presented in the left panel of 
FigureH The FIRST image (beam size = 5".40 x 5".40) 
of the background QSO J122106. 87+454852.1 with its 
lobes is shown in blue in (enlarged version in the right 
panel) and its D-configuration continuum image (beam 
size = 52". 17 x 38". 11) is shown with white contours. 
The sightlines (sightline (SL)# 12, 13, and 14) to the 
background continuum source probes the galaxy at im- 
pact parameter of 3.2, 2.8, and 2.6 kpc. The H I distri- 
bution in this galaxy is extended, extending 2-3 times be- 
yond the optical size of the galaxy. The H I is asymmet- 
ric and extends 5'. 6 to the southeast of the galaxy. The 
Westerbork observations of neutral Hydro gen in Irregu- 
lar a nd SPiral galaxies (WHISP Survey; ISwaters et al.1 
|2002[) also imaged this galaxy. The subsequent analysis 
by ISwaters et al.l (|2009t ) found that its H I distribution 
is clumpy and lacking ordered rotation. 

As we observed UGC 7408 with both the VLA and 
GBT, we were able to compare the HI emission from 
the galaxy that we obtained from each telescope. Fig- 
ure |4] shows a good match between the spectra from the 
two telescopes. The VLA recovered integrated flux of 
7.4 Jy km s _1 (removing negative pixels) as compared to 
8.1 Jy km s" 1 of GBT. This confirms that the VLA data 



recovered more than 91% of the H I flux density detected 
by the GBT. The H I spectral profile shows a deviation 
from a Gaussian profile at the peak of the emission and 
in the high velocity wing between 480 and 500 km s" 1 
. The VLA spectrum obtained from a region of equal 
size as the synthesized beam in the uniform weighting 
case (the black rectangle in the left panel of Figure [3]) 
over the quasar reveals two emission peaks with a dip at 
442.2 km s _1 . This could be a sign of absorption and is 
discussed in detail in section l4~2l 

In the other two galaxies with H I emission 
(J141629.25+372120.4 and J160659. 13+271642.6), the 
expected H I diamete r estimated using th e relationship 
shown in Figure 7 of (jSwaters et al.ll2002T ) are 30.7 and 
60.2 kpc, respectively. Thus, the sightlines, SL# 22 and 
23 are expected to pierce through the H I distribution 
in these systems. Although we detect dips in the GBT 
spectra, it is unclear if there is real absorption or are 
overlapping individual emission features mimicking ab- 
sorption. Hence, we consider these two cases as nondc- 
tections. A detailed discussion of possible beam effects 
and how the GBT spectra in these two cases may be a 
result of the beam being filled-in by emission is presented 
in Section 14.21 

Besides the H I emitting galaxies, we detected one 
unambiguous H I absorption features at 21 cm in 
the foreground galaxy GQ104 2+0747 (SL#10) a re - 
sult we discussed in detail in IBorthakur et al.l (|2010( l. 
GQ1042+0747is a dwarf spiral galaxy of luminosity 
0.048 L». The sightline pierces through the optical disk 
in this moderately star-forming galaxy. We estimated the 
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Fig. 3. — LEFT: VLA D-configuration H I 21 cm column density contour image of UGC 7408 overlaid as green contours on an SDSS-r 
band image. The contours levels arc at 1.8, 3.6, 7.2. 14.4, 28.8, &: 57.6 X 10 19 x cm -2 respectively. The beam size associated with the 
H I data obtained using natural weighting is shown at the bottom left corner (67". 51 X 45". 31). The blue contours represent the VLA 
FIRST survey image of the background quasar J122106. 87+454852.1 at 1, 5, and 25 mjy/beam and the white contours represent the VLA 
D-configuration 21 cm continuum image obtained using uniform weighting at 5 and 50 mjy/beam (Position angle= —86.75°). The higher 
resolution FIRST survey images shows three peaks from the background source, which are blended together into a single source in the 
D-configuration image. The black square represents the region corresponding to the red spectrum in Figure [4] which is roughly the size of 
the beam for the uniformly weighted data. RIGHT: A magnified version of the left hand panel showing the background quasar as imaged 
by the FIRST survey. The brightest background source at 1.4 GHz is the southern component of the quasar. The three peaks of the 
background quasar mark the sightlincs 12, 13, and 14. 



kinetic temperature of the cloud to be < 283 K based on 
the absorption line FWHM of 3.6 km s _1 . Assuming the 
spin temperature to be same as the kinetic temperature, 
we estimated the column density of the absorber to be 
JV(H I) < 9.6 x 10 19 cm" 2 . A detailed analysis of the 
properties of the absorber, the foreground galaxy, an d 
its environment is presented in iBorthakur et all ( 2010f ). 
We did not detect any H I absorber in our GBT data 
from the three H I emitting regions. Part of the problem 
may be the large GBT beam size, which can include a 
substantial portion of H I emitting regions that may fill 
up the absorption. This possibility is also discussed in 
detail in section 14.21 Table [3] presents the 3 a limiting 
optical depths for each of the sightlincs along with the 

=100 K 
1 and a line-width of 



flux density sensitivity achieved. Assuming T s 
for the absorbing gas with f c , 



■>5 km s , the corresponding limiting column densi- 



ties would range from 0.09 
median of 1.46 x 



5.83 x 10 20 cm" 2 , with a 



4. DISCUSSION 



4.1. 



Covering Fraction 

The covering fraction of cold gas in galactic halos has 
been a matter of debate. While low-ionization gas in 
the Lyman a and Mg II transition have been detected 
in the halos of various galaxies, HI 21 cm transi- 
tions have not been e xplored extensively. For instance, 
lLanzetta et al.l (|1995l ) found that a significant portion of 
the Lyman a absorbers (1/3 to 2/3) were associated with 
normal, luminous galaxies with Lyman a absorbing halos 
of size 160 /i _1 kpc. They found an anti correlation be- 
tween the strength of the absorption (equivalent width) 
and galaxy impact pa rameter, which was later confirmed 
by ot her authors (e.g. lTripp et alJl998tlChen et al.lll998l 
120011 but see also Bowen et al. 2002; Wakker & Savage 



2009). In this study, our aim is to estimate the frequency 
of occurrence of 21 cm absorbers in the halos and ex- 
tended disks of galaxies. Figure [5] shows our findings in 
terms of two plots. The left plot shows the variation 
of measured optical depth as a function of impact pa- 
rameter and the plot on the right shows the same quan- 
tity as a function of the ratio of impact parameter over 
the r-band Petrosian radius of the foreground galaxies. 
While the left plot shows the strength of the absorption 
or limiting strength versus the rest frame projected dis- 
tance, the plot on the right shows the same quantity as 
a function of the optical size of the foreground galaxy. 
In addition, the optical colors of the foreground galax- 
ies are denoted using respective colored symbols, and 
the luminosity with symbol-size, as described in the leg- 
end. Our single 21 cm absorber detection towards, SDSS 
J104257.58+074850.5, is shown as a filled circle while all 
the other non-detections are shown with open circles with 
arrows pointing down. W e do not include the l imiting op- 
tical depth values for the iGupta et al.l (|2010[ ) sample as 
their channel width is more th an an order o f magn itude 
higher than ours. In addition, IGupta et al.l (|2010f ) pre- 
sented their limiting optical depth for equivalent velocity 
resolution of 10 km s _1 . In contrast, we present our op- 
tical depth limits for a resolution of 1 km s _1 . The high 
resolution guarantees that such a study will detect ab- 
sorbers similar to the one in GQ1042+0747 or for H I at 
a temperature of ~50 K. 

The nondetection of cold H I can be understood in 
terms of three possibilities: (1) the nonexistence of cold 
gas along the QSO line of sight, (2) a low covering frac- 
tion of the background QSO by the absorber or high spin 
temperature of the absorber, and/or (3) confusing effects 
due to the large beam size of the GBT (FWHM of 9'.1). 
The possible existence of cold gas clouds with Damped 
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Fig. 4. — LEFT: A comparison plot showing the GBT HI spectrum from a region of size 9.1' (equivalent to the GBT beam at 21 cm) 
toward UGC 7408 in black and the VLA spectrum, obtained using natural weighting, encompassing the entire HI associated with UGC 7408 
in blue. The good match between the two spectra proves that the VLA recovered almost 100% of the H I associated with this galaxy. 
The VLA H I spectrum from the region marked by the black rectangle (Figure [3] right panel) of size 40"x50" is shown in red in the 
comparison plot. The spectrum was obtained from data reduced using uniform weighting, which resulted in a better spatial resolution. 
RIGHT: A magnified version of the red spectrum showing two emission peaks and a dip in between at 442.2 km s _1 . Interestingly, the 
stronger peak at 436.3 km s — 1 is toward the low-velocity end of the profile unlike the global H I profile of the galaxy at that velocity. The 
green solid line and the blue dotted line represent Gaussian fits to the emission and the absorption profile respectively. The absorption is 
narrow with FWHM of 4.75 km s" 1 corresponding to a kinetic temperature of < 493 K. 

Lyman a (DLA) column density of 2 x 10 20 cm -2 in the luminosities in our sample, n s ~2x 10 -2 Mpc -3 . Based 

on the Equation^ the dN/dz statistics require f^ A = 



region surveyed can be estimated using the cross section 
impli ed by the dN/dz of DLA systems (see e.g., section 
6.1 in lSchave et al.ll2007j ). where the covering fraction of 
such clouds is described as 



./; 



DLA 



« >< " (f ) 



10- 2 Mpc 



-3 



pc/ 



(2) 

where dN/dz is the rate of incidence of DLA sys- 
tems at z = 0, n g is the comoving number den- 
sity of galaxies affiliated with the absorbers, and r 
is the radius of the absorber cross section. dN/dz 
for 21 cm column densities comparable to those of 
DLAs in the nearby universe was computed by vari- 
ous authors from large H I surveys like the Arecibo 
Dual-Beam Surve y (ADBS) and HI Pa r kes A ll Sky 
Survey (HIPASSh [Rosenberg fc Schneiderl (120031) found 
dN/dz = 5.3 (± 1 . 3) x 10~ 2 from the ADBS w hereas 
iRvan- Weber et all (|2003l) and lZwaan eTail (J2005I ) found 
dN/dz = 5.8 (± 0.6) x 10" 2 and 4.5 (± 0.6) x 10~ 2 
respectively from the HIPASS. For our calculation, we 
adopt dN/dz = 5.3 x 10~ 2 from, Ro senberg fc Schneiderl 
which is also the median of the three values cited above. 
This value is consistent with the DLA incidence per unit 
redshift n DLA (z) = (0.044 ± 005) x (1 + z )i -27±o.ii 
for < z < 5 as described by iRao et all (120061 ) based 
on their HST survey of SDSS quasars with intervening 
Mg II-Fe II absorption systems. The number density 
of galaxies within a magnitude range can be estimate d 
from the Schechter luminosity function (|Schechterlll976l ). 
For our analysis, we use the Schechter function fit of 
</>» = 1.49 x 10~ 2 h 3 Mpc" 3 , M* - 51oginfe = -20 .44, 
and a = —1.05, estimated by iBlanton et al.l (|2003[ ) for 
SDSS r-band data k-corrected to z=0.1. For the range of 



for cross section of 14 kpc or less. The covering fraction 
drops as the cross section of the galaxy increases. For 
our sample, at a cross section radius of 55 kpc, f^ A , 
drops to «0.08, and at 90 kpc it further drops to «0.02. 
The presence of detectable 21 cm absorption not only 
depends on f^ A and the column density of the gas but 
also on the spin temperature, T s , as well as the fraction 
of the background quasar covered by the absorber. This 
fraction, which we will refer to as / hereafter, which is 
different from f^ A - For an H I absorption feature 
the column de nsity can be estimated u sing the following 
equation from lRohlfs fc Wilson I (J200Q[ ): 



A(H I) = 1.823 x 10 18 -y J T 21 (v)d% 



(3) 



where T2i(i>) is the 21 cm optical depth as a function of 
velocity in km s . In cold clouds collisional excitation is 
the primary source of excitation for the hydrogen atom 
emitting the 21 cm line. Therefore, the above equation 
can be simplified by assume the spin temperature to be 
same as the kinetic temperature, T&, which in turn can 
be related to the full-width at half maximum (FWHM) 
of the absorption feature by 

T k < 21.855 (A v) 2 (4) 

where Av is the FWHM velocity in km s" 1 . Therefore, 
Equation [3] can be rewritten as 

n3/2 



7V(H I) = 0.390 x 10 



IS 



t: 



f 



' avg 



cm 



(5) 



where r, 



lg represents the average optical depth such that 
A v = J T-2i(v) dv . Figure IH1 shows the vari- 
ation of H I column density as a function of optical 



i avg 
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Fig. 5. — LEFT: Impact parameter, p, versus the GBT integrated optical depth over channels corresponding to 1 km s" 1 (i.e. 3 channels). 

J v P " ak —0 5 T d v > f° r our sample. The filled circles indicate 21 cm absorbers and open circles represents the 3 a limiting values. The colors 

show color classification of the galaxy as discussed in Section 12.21 The peak optical depth per km s — 1 for the absorber detected in the 
VLA data of UGC 7408 is shown as a filled star. RIGHT: Ratio of impact parameter to optical Petrosian r-band radius, p/R op t, versus 
the integrated optical depth. The dotted line represents p/R pt = 1. The sightlinc probing GQ1042+0747 is the only sightline to probe 
the optical disk. This is also our only unambiguous 21 cm absorption detection with the GBT. 



depth, assuming / — 1. The contours represent log- 
arithms of H I column density and are labeled accord- 
ingly. A column density of 2 x 10 20 cm -2 , corresponding 
to a DLA system, is shown with a solid line with hllcd 
circles overplotted and labeled 20.3. The filled cyan cir- 
cle shows the absorber associated with GQ1042+0747. 
The dashed line represents the ma ximum spin temper- 
ature measured in the Milky Way (jBraun fe Walterbosl 
Il992t ). although the majority of their observations for the 
Milky Way and M31 have T s below 300 K. The absorber 
associated with GQ1042+0747 is consistent with these 
observations. However, higher values of spin tempera- 
ture have been commonly seen at higher redshifts and 
even in lo w-redshift absorbers associated with non-spiral 
galaxies (lLane et al.l 119981: IChengalur fc Kanekarl 120001: 
iKanekar fc Chengalurl 12001. 200jJ. Our survey aims to 
detect cooler clouds with DLA equivalent column densi- 
ties up to 1100 K (if / = 1). It is clear from Figurc[6]that 
even for our best limiting optical depth value of 0.0139, 
if T s > 1100 K then the corresponding H I column 
density would have to be well over the DLA limit to be 
detectable as a 21 cm absorber. On the other hand, we 
can be certain that a non-detection confirms that there 
are no cold clouds with T s < 100 K (assuming / = 1) 
and column density equivalent to a DLA in our entire 
sample, with the exception of our noisiest sightlinc (SL^ 
18 With Tavg < 0.64). 

4.2. Cold gas in the Dwarf Galaxy UGC 7408 

The GBT spectrum towards 122106.854+454852.16 
does not reveal any absorption from UGC 7408 despite 
the fact that there is H I emission in the region in front 
of the quasar as seen in the VLA D-configuration im- 
age (see Figure |3]). The H I column density in the 
foreground region within the two lobes is between 1.44 
and 2.88 x 10 20 cm" 2 at the resolution of VLA D- 
configuration. In this section, we investigate the possible 
causes that may be responsible for the lack of absorp- 
tion features in the GBT spectrum and what that may 
mean for the rest of our sample. The nondetection of cold 



gas in absorption, despite the good sensitivity achieved, 
could be due to any of the three following underlying 
causes. For an H I emitting galaxy such as UGC 7408, 
it is highly possible that the absorption is filled in by 
emission from the rest of the galaxy. This also extends 
to the other two H I emitting foreground galaxies in our 
sample. The second possibility is that the H I is mostly 
warm and does not absorb the background continuum. 
The third possibility is that the ISM in this galaxy is 
highly patchy and N(H I) is much lower along the QSO 
sightlines than the average inferred at larger scales from 
the emission map. 

Ideally, 21 cm absorption studies require high angu- 
lar resolution to prevent absorption from being filled in 
by em ission from the rest of the system. iDickev et al.l 
(|2000T l suggests a resolution of < 10" for extragalactic 
sources. Unfortunately, the high sensitivity and spec- 
tral resolution achieved by the GBT comes at the cost 
of spatial resolution. On the other hand, our VLA D- 
configuration beam is much smaller than the GBT beam 
and this provides us with the opportunity to search for 
absorption from the region on the line of sight to the 
background quasar. The right panel of Figure |4] shows 
the VLA spectrum (in red) from the region marked by 
the black rectangle in Figure [3l which is also the size 
of beam for the uniform weighted data. The spectrum 
exhibits peaks at 436.3 and 449.3 km s" 1 with a dip in 
between at 442.2 km s _1 . The fact that the absorption 
dip sits right at the peak of the local emission suggests 
that the dip is an absorption feature. Furthermore, the 
ratio of the strength of the two peaks in the red spec- 
trum is opposite to that of the global H I profile of 
UGC 7408 at that velocity. This further suggests that we 
are spatially and kinematically resolving the underlying 
H I emission. Even with the VLA D-configuration beam, 
we are still primarily detecting the emission from the 
foreground region. For an H I cloud and a background 
quasar contained within a single beam, the observed line 
temperature, Tl, can be expressed as Equation (12.22) 
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of iRohlfs fc Wilson I (t2000h 

T L = (f cl T s - f f c T c )(l - e- T ), (6) 

where, f c i and /c are beam filling factors of the H I cloud 
and continuum source, respectively, /o is the fraction 
of the continuum source covered by the H I cloud, 
and Tc is the brightness temperature of the background 
source. From the FIRST image of background quasar, 
J122106.854+454852.16, we find that f c < 0.33 and 
Tc = 25 K for the VLA in D-configuration, and thus 
we expect the emission to be significantly stronger than 
absorption even for very cold clouds (T s ~ 25 K). 

In order to get a proper estimate of the H I column 
density associated with this dip, it is necessary to sepa- 
rate the emission profile from the absorption. By assum- 
ing the emission from the region covered by the VLA 
beam is a Gaussian, we fit the emission profile (exclud- 
ing the dip). The green solid line in the right panel of 
Figure |4] shows the emission profile. Similarly the absorp- 
tion profile is also fitted with a Gaussian profile yield- 
ing a FWHM=4.75 km s" 1 centered at 442.2 km s" 1 
. The column density associated with this profile can 
be estimated using equation [3] If we integrate this pro- 
file and divide by the total background flux density of 
67 mJy (background continuum flux density of 56 mJy 
+ background H I emission flux density of 11 mJy), 
we get N(HI)= 3.8 x 10 19 (10 ^ // cur 2 . The ki- 
netic temperature associated with this profile of FWHM 
=4.75 km s" 1 is < 493 K and if T s = T kl then the col- 
umn density associated with this absorber can be as high 
as 1.86 x 10 20 cm" 2 (assuming / = 1). The VLA D- 
configuration data may also suffer from emission filling-in 
the absorption as the beam encompasses a much larger 
region than the background continuum peaks. Higher 
resolution VLA A-configuration and/or VLBA data are 
needed to resolve the structures in the background quasar 
and to map the H I only at the region in front of the 
continuum source. These observations will be crucial in 
understanding the small-scale structure of the ISM of 
UGC 7408. 

In the absence of higher resolution data, it cannot be 
proven that two peaks and the dip in the VLA spectra 
represents a signature of absorption. The feature may 
be part of the emission profile itself. If the feature were 
not an absorber, then it would mean that the HI is at 
a higher temperature than that seen in the Cold Neutral 
Medium (CNM) of the Milky Way. Spin temperatures of 
around 1000 K or higher have bee n seen at low-redshifts 
(z < .5) in non-sp i ral sy stems (|Kanekar fc Chengalurl 
12001 . iDickev et all (f2000h found that in the SMC the 
fraction of cool-phase H I is much lower than the overall 
H I abundance. The ISM properties of dwarf galaxies 
arc different from bigger galaxies like the Milky Way and 
this may be the reason for the low cold-H I abundances 
in their ISM. Dwarf galaxies have lower abundances of 
metals and thus, cooling from fine-structure lines (espe- 
cially from carbon and oxygen) in these galaxies is much 
less predominant than brighter galaxies (see Dickey et 
al. 2000 for details). In addition, the typical dust to 
gas ratios in dwarf galaxies are much lower than in the 
Milky Way and thus dwarf galaxies are inefficient in con- 
fining the UV-ionizing radiation to individual H II re- 
gions. Therefore, UV radiation can effectively heat up 
the cooler gas in dwarf galaxies to a larger radii outside 
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Fig. 6. — Variation of H I column density as a function of average 
optical depth, T avg and spin temperature, T s for absorbers with 
covering fraction, / = 1. The contour levels mark column density 
in logarithmic scale. The contour at log N(HI)= 20.3, which is 
indicated by a solid line with filled circles, represents the DLA limit 
of N(HI)=2xl0 20 cm -2 . The horizontal dashed line represents the 
maximum T s detected in our Ga laxy within the range of t values 
shown in the plot (see Fig 2a of Braun & Waltcrbos 1992). The 
maximum T 3 for M31 is around 1000 K, although the majority of 
the data for both M31 and the Milky Way lie below 300 K. Our 
21 cm absorber detected in the foreground galaxy GQ1042+0747, 
is plotted with a cyan circle. This is consistent with absorbers 
detected in the Milky Way. The GBT 3<r limiting optical depths for 
the rest of our foreground galaxies that show no 21 cm absorption 
are plotted at T s = 300 K with arrow headed triangles and span a 
wide range of column densities. 

the stellar disks. Hence, due to increased heating and 
reduced cooling, it is very likely that most of the H I in 
dwarf galaxies is at a much higher temperature. 

Another possibility for the nondetection of H I ab- 
sorption is the lack of H I in the foreground region. This 
may be a result of a clumpy ISM, which when observed 
with VLA D-confi guration resolution get smoothed out. 
iPuche et al.l (|1992T ) found holes in the H I distribution of 
Holmberg II, a much-studied star-forming dwarf galaxy, 
with diamete r s rang ing from ~ 100 to 1700 pc. Similarly, 
iBegum et al.l (|2008l) also found gaps and/or holes in the 
H I distribution for dwarf galaxies such as DD043. Maps 
from the WHISP survey also show that galaxies have 
similar patches of low intensity H I emission even in the 
inner regions. In the case of UGC 7408, since we have 
three sightlincs through the ISM at impact parameters 
2.6, 2.8, and 3.3 kpc, according to the last possibility, the 
ISM needs to be extremely clumpy. Also it requires the 
clumps to have very high column densities so that the av- 
erage column density in emission when smoothed to VLA 
D-configuration resolution is around 2 x 10 20 cm~ 2 . In 
conclusion, we emphasize the need for high spatial resolu- 
tion VLA A-configuration imaging. These observations 
will not only provide information on the dumpiness of 
the ISM in dwarf galaxies but also will shed light on 
the spin temperature of the H I in such systems. Fur- 
thermore, the background QSO is a UV bright source 
and can be observed with the Cosmic Origin Spectro- 
graph (COS) on the Hubble Space Telescope (HST). The 
spectroscopic UV data would provide complementary in- 
formation which will allow us to learn more about the 
multiphase ISM including the warmer phases of neutral 
gas in UGC 7408. 
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4.3. Connection of 21 cm Absorbers to ISM of Galaxies 

Our only unambiguous GBT detection comes from the 
foreground galaxy GQ1042+0747, which is intercepted 
by a QSO sightlinc that has the smallest impact parame- 
ter of all the QSO-galaxy pairs in our sample, and which 
pierces the stellar disk of the galaxy. As discussed in 
iBorthakur et al.l (|2010t) , the absorber is very narrow and 
has a column density of iV(H I) < 9.6 x 10 19 cm~ 2 , 
i.e., lower than a DLA. It is evident from our sample 
that we do not find cold clouds with column densities 
of the order of DLAs or more, associated with the halos 
of galaxies ubiquitously. The second possible absorber, 
associated with UGC 7408, probably arises in the bound 
ISM of the dwarf galaxy. 

Other studies have found 21 cm absorbers and some of 
them have b een associated with tidally stripped matter. 
For example, IStocke et al.l (|1991l ) found a 21 cm H I ab- 
sorber associated with a tidal feature of NGC 3067. Cool 
gas associated with galactic outflows has been observed 
in the h alos of galaxies probed by various metal-line tran- 
sitions. IVeilleux et al.l (J2005I ) observed such outflows in 
nearby star burst galaxies a nd iBland-Hawthorn fc Cohen! 
(|2003f ) and iKeenev et al.l (|2006t) found indications of a 
bound outflow in the central region of the Milky Way. 
However, the nondetection of cold gas in the halos of 
galaxies in our sample suggests that cither the amounts of 
H I in these systems are not high enough or T s > 1100 for 
the outflowing material. Since most of the measurements 
of outflowing material come from metal-line transitions, 
the temperature of the gas may be well above our limit 
to detect 21 cm H I absorption even for our best-case 
sensitivity scenario. In addition, the radio background 
quasars are often extended and thus a clumpy structure 
of out-flowing material can have a relatively low cover- 
ing fraction, /, thus significantly increasing the column 
density associated with our limiting optical depth. 

The same argument can also be extended to accret- 
ing material or cold clouds condensing out of the halos. 
The above conclusions hold irrespective of the metallic- 
ity of the cold clouds, and therefore provide information 
about accreting and outflowing material with the same 
accuracy. Out of 15 foreground galaxies in our sam- 
ple, three of them have masses < 10 103 M© with no 
detectable H I emission, and sightlines toward the back- 
ground quasars through the halo of the galaxies. How- 
ever we did not find any 21 cm absorber in their halos. 
This is not surprising as "cold" accretion could easily 
have spin temperatures that makes it difficult to detect 
in 21 cm absorption. We also do not find ubiquitous 
21 cm absorption in the halos of luminous galaxies un- 
like Lyman a studies. 

5. SUMMARY 

The process of gas flow into galaxies and subsequent 
condensation into cold gas has been a topic of various 
observational and numerical studies. In this paper, we 
investigated the presence of cold gas in the halos and 
extended disks of low-redshift galaxies. We find the oc- 
currence of H I 21 cm absorption is rare (1/12) as we 
move outwards from the optical disk of the galaxies for 
our sample. We exclude the three, H I-emitting fore- 
ground galaxies as nondetections as any absorption along 
the background lines of sight may have been filled in by 
emission from H I in the ISM of the galaxy itself. The 



only case where we unambiguously discovered atomic gas 
in absorption was in the stellar disk of the dwarf galaxy 
GQ1042+0747 at an impact parameter of 1.7 kpc from 
the QSO sightline SDSS J104257.58+074850.5. For sys- 
tems with p < 7.5 kpc, which was reported as the me- 
dian impact parameter for H I systems corres ponding 
to DLA column densities bv lZwaan et al.l (|2005l ). we de- 
tected H I in one out of two galaxies. Our detection 
rates for sig htlines with p < 15 kpc do not match the 
results from IGupta et al.l (|2010D . This may be partly 
because of the inclusion of pairs selected due to the 
presence of Ca II and Na I absorp tion (sightlines from 
iCarilli fc van Gorkomlll99l in the IGupta et~afl ff20T?)h 
sample. 

On the other hand, we detected HI in emission 
in three foreground galaxies including the dwarf galaxy 
UGC 7408. The H I masses range from 0.9 to 
81.3 x 10 8 M . The VLA H I imaging of UGC 7408 
shows that its H I envelope extends to 2-3 times the opti- 
cal size of the galaxy in the form of an ellipsoid extending 
in the southeast direction. Interestingly, the background 
continuum sources (SL# 14, 13, and 12) at impact pa- 
rameters of 2.6, 2.8, and 3.3 kpc from the optical center 
of the galaxy is well within the H I distribution seen in 
VLA imaging. We did find a dip in the H I emission at 
the position of the background quasar. The dip is cen- 
tered at 442.2 km s" 1 and has a FWHM of 4.75 km s" 1 
corresponding to a kinetic temperature of < 493 K. As- 
suming the spin temperature to be same as the kinetic 
temperature, we estimate the column density associated 
with this dip to be 1.86 xlO 20 cm~ 2 . However, higher 
resolution imaging is required to confirm this detection, 
because the dip may be a result of the emission profile 
from the foreground ISM of UGC 7408 and not an ab- 
sorption signature. This may suggest that the H I seen 
in the VLA imaging may have too high a temperature to 
detect 21 cm absorption or the H I distribution in this 
galaxy is highly patchy with no H I structures at the 
position of the background continuum source. 

Our study was designed to be most sensitive to cold 
gas with temperatures observed in the CNM of the 
Milky Way. However in a few cases our sensitivities 
are high enough to look for DLAs at temperatures 
above 500 K. For instance, along the sightline where we 
achieved the best 3 a limiting optical depth of 0.01, we 
were sensitive to DLA systems even at >1000 K. On the 
other hand, for low temperature systems, this means 
that we could detect absorbers with column densities as 
low as 5.6 x 10 18 cm" 2 if T s = 100 K. Therefore, we 
can conclude that in most cases, there is no H I with 
DLA column densities at temperatures similar to that of 
CNM in Milky Way in the halos and extended disks of 
galaxies. Although our data are insufficient to confirm 
the inflow and/or outflow of neutral gas through the 
galactic halos, our nondetection of H I in the halos 
does hint that the process of condensation of warmer 
gas into H I may be occurring at regions closer to the 
optical disk of galaxies. Our present understanding of 
the cold gas distribution in the halos and extended disks 
of galaxies is quite limited. To extend our knowledge, 
high-resolution H I imaging of galaxy-quasar pairs 
as well as complementary information about warmer 
phases of gas through UV observations is crucial. In 
the future, the Expanded VLA (EVLA) with broad 
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bandwidth and the COS with high sensitivity will be 
able to probe systems through observations that would 
yield information on the physical nature and distribution 
of various gaseous species in the halos and extended 
disks of galaxies. 
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